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ff^ ■ mixing, while the standard model prediction is very small. Such a large phase may 

O . 

fSJ , imply sizable new physics effects in i?g — B^ mixing. We compute the gluino-mediated 

' supersymmetry contributions to B^ — B^ mixing, Bg — ?• and B — t- Xg^ 

decays in the frame of the mass insertion approximation. Combining the constraints of 
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^ . constrained LL and RR insertions in Bs — )• K^*)~ K^*")^ decays are small because of the 
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. . .' absence of gluino mass enhancement. For ??T,|/m? = 9, the constrained LR insertion 

can provide sizable contributions to all observables of Bg — )• decays except 

A'q^p{Bs — 7- K~ K^)^ and many observables are sensitive to the modulus and the phase of 
the LR insertion parameter. Near future experiments at Fermilab Tevatron and CERN 
LHC-b can test these predictions and shrink/reveal the mass insertion parameter spaces. 
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1 Introduction 



The flavor changing neutral current (FCNC) processes in the b s transition are sensitive to 
the effects of New Physics (NP) beyond the Standard Model (SM). Recently, both CDF and 
D0 collaborations have announced their measurements of extracted CP violating phase (pl^'^^ 
associated with — B'^ mixing 1 j3||. The CP violating phase measured by both CDF and D0 
is (pi^^'^ E [0.20, 2.84] at 95% C.L. J, which is much larger than its SM value 0W'SM ^ 2/3™ = 



2arg 



VtsV, 



tb. 



0.04 



a 



More recently, the D0 collaboration has reported evidence for an 



anomalously large CP violation in the like-sign dimuon charge asymmetry in semileptonic B- 
hadron decays = (-9.57 ± 2.51 (stat) ± 1.46(syst)) x 10"^ (loL which differs by 3.2 standard 
deviations from the SM prediction A^f^ = (-2.3t[JJ) x 10"^ p, |ll|. Although the errors of 
the data are still large, these deviations from the SM could be attributed to the presence of 
non-SM flavor violationin the b ^ s,d nonleptonic decays. 

Recently, the CDF collaboration has made the first measurement of charmless two-body 



m 



13|. The measurement 



B, K~K+ decay, ^ K-K+) = (26.5 ± 4.4) x lO^*' 

is important for understanding Bg physics, and also implies that many Bs decay modes could 
be precisely measured at the LHC-b. Comparing with the theoretical predictions for these 
observables in Refs. 14Nl6|. one would find that the experimental measurements of branching 



ratio are in agreement with the SM predictions within their large theoretical uncertainties. 
However, NP effects would be still possible to render other observable deviated from the SM 
expectation with the branching ratios nearly unaltered |17| . 

Supersymmetry (SUSY) is an extension of the SM, which emerges as one of the promis- 
ing candidates for NP beyond the SM. In general SUSY, a new source of flavor violation is 
introduced by the squark mass matrices, which usually can not be diagonalized on the same 
basis as the quark mass matrices. This means gluinos (and other gauginos) will have flavor- 
changing couplings to quarks and squarks, which implies the FCNCs could be mediated by 
gluinos and thus have strong interaction strength. It is customary to rotate the effects so they 
occur in squark propagators rather than in couplings, and to parameterize them in terms of 
dimensionless mass insertion (MI) parameters {S^B)ij with {A,B) = {L,R) and (i,j = 1,2,3). 

-B^ — -Bg mixing, Bg — t- and B — )■ Xs'j decays are all induced by the b ^ s transi- 

tion, and they involve the same set of the MI parameters. Inspired by the recent measurements 
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from CDF and D0 collaborations , we study B2 — -B° mixing, Bg — K^*^ and B — 



decays in the usual MI approximation (l8l . |19| of general SUSY models, where flavor violation 
due to the gluino mediation can be important. The chargino-stop and the charged Higgs-top 
loop contributions are parametrically suppressed relative to the gluino contributions, and thus 



are ignored following 



19 



3I-I22I ■ Following the similar way to our previous article 23|, we consider 
the LL, RR, LR and RL four kinds of the Mis with m?/m? = 0.25, 1, 4, 9, respectively. We find 
that the LL and RR insertions for all cases of m?/m? values as well as the LR insertion for 
m?/m? = 9 case could explain current experimental data simultaneously. For m?/m? = 9, the 
constrained LR MI could significantly affect all observables of Bs — J- K^*^~ K^*^~^ decays except 
Aq<p{Bs — )■ K~K'^) without conflict with all related data at 95% CL. While the constrained 
LL and RR insertions from B^ — B^ mixing have small effects in Bs — ?■ K^*^~ K^*^'^ decays 
because of the absence of gluino mass enhancement. Therefore, with the ongoing 5-physics at 
Tevatron, in particular with the onset of the LHC-b experiment, we expect a wealth of Bs data 
and measurements of these observables could restrict or reveal the parameter spaces of the LR 
(LL and RR) insertions in the near future. 

The paper is arranged as follows. In Sec. 2, the relevant formulas for Bs — )■ K^*"^' K^*"^^ 
decays and 5° — 5° mixing are presented. Sec. 3 deals with the numerical results. Using 
our constrained MI parameter spaces from Bs -> decay and -B^ — 5° mixing, we 

explore the MI effects on the other observable observables, which have not been measured yet 
in Bs — )■ K^*'^~ K^*'^^ decays. Sec. 4 contains our summary and conclusion. Theoretical input 
parameters are summarized in the Appendix. 

2 The theoretical frame 

2.1 B, K^*yK^*)+ decays 
2.1.1 The decay amplitudes in the SM 

In the SM, the effective Hamiltonian for the b — ?■ suu transition at the scale /i ~ rrih is given 



by M 

n'JjiAB = 1) = % E f c^f ^g? + ct'Ql + E crQ. + c^i'Qi, + cl^Q.^ +h.c., (i) 
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where \p = VpbV* with p G {u, c} are Cabibbo-Kobayashi- 



coefficients within the SM Cf *^ can be found in Ref. 
given as 



Vlaskawa (CKM) factors, the Wilson 



24| . and the relevant operators Qi are 







Q^ = 


{Pal''Lb[s){s,37i,Lp^), 


Q3 = 


q' 


Q4 = 


{sp-f^Lbo,)J2iqali^Lq'p), 

q' 


Q5 = 


q' 


Q6 = 


{s^YLb.)Yiqai,H)^ 

q' 


Q7 = 


3 

q' 


Q8 = 


3 

3' 


Q9 = 


3 

-{Sal'^Lba) eg' {q'pll.Lq'13) , 


Qio = 


3 

-{SfS-f'^Lbc,) Y (^q'iicnMfi 

q' 


Qr-y = 


g 




ll^m,-s^a^''RT:^b^G%, 


and /3 are color indices, and L{R) 


= (1T75). 



(2) 



With the effective Hamiltonian given in Eq. ([T]), one can write the decay amplitudes for the 
relevant two-body hadronic B — ?■ M1M2 decays as 



A^^{B ^ MiM2) = [MiM2\'Hlfj{^B = l)\B] 



(3) 



The essential theoretical difficulty for obtaining the decay amplitude arises from the evaluation 
of hadronic matrix elements {MiM2\Qi{^)\B) , for which we will employ the QCD factorization 
(QCDF) 25| throughout this pap er. We will use the QCDF amplitudes of these decays derived 



in the comprehensive papers [15 



26| as inputs for the SM expectations. 



2.1.2 SUSY effects in the decays 

In the SUSY extension of the SM with conserved R-parity, the potentially most important 
contributions to Wilson coefficients of penguin operators in the effective Hamiltonian arise 
from strong penguin and box diagrams with gluino-squark loops. They contribute to the 
FCNC processes because the gluinos have flavor- changing coupling to the quark and squark 
eigenstates. In general SUSY, we only consider these potentially large gluino box and penguin 
contributions and neglect a multitude of other diagrams, which are parametrically suppressed 
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Figure 1: Penguin diagrams for b suu process with gluino exchanges at the first order in 
mass insertion, where h,k,m = L,R. 



(a) 



















1 / 












(b) 



Figure 2: Box diagrams for b — t- suu process with gluino exchanges at the first order in mass 
insertion, where h,k,m = L,R. 



by small electroweak gauge coupling 19l-|22||. The relevant Wilson coefficients of the b suu 
process due to the gluino penguin or box diagrams, which are shown in Fig. [1] and Fig. |2l 
respectively, involving the LL and LR insertions are given (at the scale /i ~ mw ~ rrig) by 



27|-|29| 



'SUSY, 



C 



^SUSY^ 



^SUSY^ 



2V2GF\tm} V 9 
'2V2GFXtml V 3 



a 



<SUSY, 



rrin] = 



23, 



2V2GFXtml V 3 



87ras{mg) 
9V2GF\tml 



(5^J23M4(X) - fc)23 (^j ^B,{x] 
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SUSY, 



ma 



+ (<^Li?)23 



m?. 



2(3; 



(4) 



where x = rn?-/rn?-^ and the loop functions Bi{x), Pi{x), Mi{x) can be found in Ref. [2 



For the RR and RL insertions, we have additional operators Q 



4=3. ..6,77,85 



that are obtained by 



L o i2 in the SM operators given in Eq. ([2]). The associated Wilson coefficients C, 



SUSY 
i=3...6,77,i 



are 



determined by the expressions as above with the replacement L ■v^ R. The remaining coefficients 
are either dominated by their SM (C1.2) or electroweak penguins (C7,,,io) and therefore small. 



The SUSY Wilson coefficients at low energy Cf ^'^^(/i ~ rrib) can be obtained from Cf ^'^^( 



ma 



in Eq. (jl]) by using the renormalization group equation as discussed in Ref. 



24| 



(5) 



where C is the 6x1 column vector of the Wilson coefficients and U^{fi, m^) 2J| is the five-flavor 
6 X 6_evolution matrix. The coefficients Cf^^'^^ and C^^^^ at the yU ~ m^ scale are given by 



30 



3l| 



riSUSY 
riSUSY 



2rySUSY , 



77 

ySUSY 



ma 



SUSY, 



m. 



(6) 



with Tj 



2.1.3 The total decay amplitudes 

For the LL and LR insertions, the NP effective operators have the same chirality with the SM 
ones, so the total decays amplitudes can be obtained from the SM ones in Refs. 
replacing 
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m by 



^SM ^ ^.M ^ ^. 



SM 



-iSUSY 



(7) 



For the RL and RR insertions, the NP effective operators have the opposite chirality with the 
SM ones, and we can get the corresponding decay amplitudes from the SM decay amplitudes 



by following replacements 32 1 



^SM ^ ^^M _ 



SM 



'-iSUSY 



(8) 



for A{Bs K-K+) and Ao,||(fis ^ K*-K*+), as well as 

CSM ^^SM^^^USY^ (9) 

for A{Bs K*"K+), A{Bs K^K*+), and A±{Bs K*-K*+). 
Then the total branching ratio reads 



B{B, ^ M1M2) = \A{Bs ^ MiM2)r , (10) 

where tb^ is the B^ lifetime, \pc\ is the center of mass momentum in the center of mass frame 
of Bs meson. 

In Bs — )■ VV decay, the two vector mesons have the same helicity, therefore three different 
polarization states are possible, one longitudinal and two transverse, and we define the cor- 
responding helicity amplitudes as ^o,±- Transverse (vAy^x) and helicity {A±) amplitudes are 
related by -4.||,_l = "^^^^ ■ Then we have 



1^(5, ^ vv)f = \Ao\' + i^+r + i^-T = i^r + i^iip + i^±r. (n) 

The longitudinal(transverse) polarization fractions fiifi) are defined by 

hABs ^VV) = ^^= „ (12) 



For the CP asymmetries (CPAs) of Bs meson decays, there is an additional complication 
due to 5° — 5° mixing. There are four cases that one encounters for neutral Bg decays, as 



discussed in Ref. 



33 



Case (i): 5^ f,B^ ^ /, where / or / is not a common final state of B^ and B^, for 
example E° ^ K~tt+ , p+ , K*~tt+ , K*^ p+ . 

Case (ii): B'^ (f = f) -^r- B^ with f-^^ = ±/, involving final states which are CP 
eigenstates, i.e., decays such as 5° — )■ K^K^. 

Case (iii): B'^ ^ (f = f) ^ 5° with f-^^ 7^ ±/, involving final states which are not 
CP eigenstates. They include decays such as 5° — )■ VV, as the VV states are not CP 
eigenstates. 
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• Case (iv): — ?■ {f&if) ^ with f^^ ^ /, i.e., both / and / are common final states 
of 5° and 5°, but they are not CP eigenstates. Decays K*-K+,K-K*+ 
belong to this case. 

For case (i) decays, there is only direct CPA (Af^p) since no mixing is involved for these 
decays. For cases (ii) and (iii), their CPAs would involve B^ — B^ mixing. The A'^'p and the 
mixing- induced CPA (A^p) are defined aJl] 

(5° ^ /) = (5° ^ /) = (13) 

\-^k\ + ^ \^k\ + ^ 



where k = 0, ||,_L for i? — )■ VV decays and A; = for B — )■ PP,PV decays, in addition, 

= CP ^^^^ = f tiS^/) CP ^^^^^ ("i) 



Case (iv) also involves mixing but requires additional formulas. Here one studies the 



our 



time-dependent decay widths for B°{t) /, fi°(t) ^ /, ^ / and / |33|-|36|. 

These time-dependent widths can be expressed by four basic matrix elements 



g = (/|?/e//|S°), h={f\n.fj\B',), 

g = h={f\neff\B''s), (14) 

which determine the decay matrix elements of B^ — > /&/ and of B^ — )■ /&/ at t = 0. We will 
also study the following observables 

/) = ^ft^, ^ /) = (15) 



^^r(5."&i?s" ^ /) = ^^r (i?°&5° ^ /) = (le) 

with A'fc = lih/g) and A'^' = lig/h). The signature of CP violation is r(i?°(t) /) ^ 
r(S°(t) -> /) and r(50(t) -> /) ^ r(50(t) -> /), which means that ^^^'"^(EOfc^o ^ /) 
-^^f (50&50 ^ /) and/or (i?o&fiO ^ /) ^ -^^^ (i?o&50 ^ /). 



2.2 5,^ - 5^ mixing 

The most general -8° — 5° mixing is described by the effective Hamiltonian 



He//(A5 = 2) = c:q: + a'q: + h.c, (i?) 



^ We use a similar sign convention to that of 37j for self-tagging and charged B decays. 



with 



Q'2 = {sPLb)i{sPLb)u 

g'3 = isPLb)sisPLb)s, 

= isPLb)^{sPJ,b),, 

Q's = {sPLb)8{sPRb)s, 



(18) 



where P, 



L{R) 



[1 =1= 75)72 and the operators Q'123 obtained from Q'123 by the exchange 



L ^ R. The hadronic matrix elements, taking into account for renormahzation effects, are 
defined as 



\B'.) 
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^Bsfk^BsB2{fJ^), 



Y^^Bjk^B^B^ifl), 
T^^Bjk^BsBiifi), 



6 



^Bsfk^BsB5{f^), 



(19) 



with 5*^^ — I — 

The Wilson coefficients receive contributions from both the SM and the SUSY loops: 
Cj' = C^^^^ + C'f^^^ . In the SM, the t — W box diagram generates only contribution to the 
operator Q[, and the corresponding Wilson coefficient C[^^^ at the nib scale is 24 1 



(20) 



where Xt = mf/m'^ and 772^ is the QCD correction. 

In general SUSY models, there are new contributions to — mixing from the gluino- 
squark box diagrams, which are shown in Fig. |3l and the corresponding Wilson coefficients 
(j/^susY ^Yie rriq scale) are given by jl9l-l22| 

„2 



Cf^^^(m,-) = --^(24x/6(x) + 66/6(x))(<^i^)i3, 



a: 



216m? 



235 



- ^- 

(a) 







1 
1 

1 " 


1 
1 

1 
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Figure 3: Feynman diagrams for — mixing in mass insertion, where h, k,l,m = L, R. 



216m|" 



216m| 



216m? 



(504x/6(x) - 72/6(x)) {SiMSU)23 - l32hix){6U)23{d1lL)23] , 
(24x/6(x) + 120/6(rE)) (51^)23(5^^)23 - l80Ux){5'iM4L)23\ (21) 



The loop functions feix) and /6(x) can be found in Ref. [27|. Other Wilson coefficients C^{'^2^3^^ 
are obtained from C'f^^"^ by exchange of L ^ R. 

The SUSY Wilson coefficients at the nib scale Cf^^^ {nib) can be obtained by 



a.(m,) = Y.T. {t'^ + Vcf'^^) r/'-C.(m,), 



(22) 



where rj' = as{niq)/as{nit). The magic number aj, fe-'^''^'* and c-^''^'' can be found in Ref. 
Renormalization group evolution of Ci,2,3 can be done in the same way as for Ci^2,3- 
In terms of the effective Hamiltonian in Eq. (ITTj) . the mixing amplitude M12 reads 

{B^,\H^„{AB = 2)\B',) 



38 



Ml 



12 



2m f 



(23) 



In the SM, the off-diagonal element of the decay matrix r^2^^ may be written as [39 1 



■^S.SM 

12 



Gpnil 



(KsK;)' fG(x,)(i?°|Qi|50) + G2(x,)(5°|g2l4°) + VT^^Ji/J , (24) 



where Xc = ml/ml, G{xc) = 0.030, and G2{x 



-0.937 at the scale 



39j. The 1/mb 



corrections Si/m are given in Ref. 40|, and l/mf corrections are not considered since they are 
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small 4l|. It is important to note that, SUSY contributions can significantly affect M{*2, but 



have little effect on which is dominated by the CKM 



hence = holds as a good approximation 

In general, the relevant CP violating phase between the — amplitude and the am- 
plitudes of the subsequent 5° and B^ decay to a common final state could be expressed as 



43|. 



avored b — )■ sec tree-level decays. 



44| 



M- 



arg 



12 



(25) 



5|. The same additional contribution 
> = + ^NP_ pg^gg sizable NP 



\ 12 / 

The SM prediction for this phase is tiny, (p^^ ^ 0.004 
due to NP would change this observed phase, i.e., 1 
contributions, the following approximation is used: 0^/^^ 

In this work, besides the CP violating phase the experimental bounds of the following 

observables will be considered: 



the Bg mass difference: AMg = 2 \M^2\ j 

3: AF, = 



the B^ width difference 



AM. 



2|Fj2|cos0,; 



the semileptonic CP asymmetry in Bs decays [46|, |47[ : Agj^ = Im 



tan05. 



3 Numerical results and analysis 

Now we are ready to present our numerical results and analysis. First, we will show our 
estimations in the SM with the theoretical input parameters listed in Table |4] of Appendix. 
Then, we will consider the SUSY effects with LL, RR, LR, and RL four kinds of the Mis and 
constrain the relevant MI parameters with the experimental data of Bs — )■ K~ , B — )■ 
and 5° — i?^ mixing. In each of the MI scenarios to be discussed, we will vary the Mis over the 
range 1(51^)231 < 1 to fully map the parameter space. We will consider the weak phases resided 
in the complex MI parameters (5^^)23 and appeared in the SUSY Wilson coefficients in Eq. (j3]) 
and Eq. (12T|) . and these weak phases are odd under a CP transformation. Using the constrained 
parameter spaces, we will give the MI SUSY predictions for the branching ratios, the CPAs 
and the polarization fractions, which have not been measured yet in Bg — j- K^*^^ K^*^'^ decays. 
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The numerical results in the SM are presented in second column of Table [H For the decays, 
the detailed error estimations corresponding to the different types of theoretical uncertainties 



have been already studied in Refs. 



in Refs. 



15 



26 



15 



26 



48 



and our SM results are consistent with the ones 



48| . For B'^ — 5° mixing, ^^Z*^"^ and Ag^ are precisely predicted in the SM, and 



the uncertainties of AMg and AF^ mainly arise from the nonperturbative quantity fssyBss 
and the CKM matrix elements. 

Table 1: The theoretical predictions for Bg — )■ decays and 5° — 5^ mixing based 

on general SUSY models with LR MI and x = 9. B and Ag^^ are in units of 10"^ and 10~^, 
respectively. The corresponding SM predictions and relevant experimental data are also listed 
for comparison. 



Ubseivabies 


Experimental ranges 
at 9070 O.L. 


biVl predictions 


SUSY values with 
i"LR)23 tor x = 9 




[17. 5a, io.UiJ 


[13.66, ^4.ozJ 


[1 ^ .53, 18.01J 


AT 


fO 05 n 3,Sl 


fO 10 211 


fn 10 211 


4>s 


[0.16,2.84] 


[0.034,0.038] 


[0.16,0.52] 




[-0.04,2.96] 


[0.02,0.05] 


[0.11,0.46] 


B{Bs -> K-K+) 


[17.70,35.30] 


[9.20,45.52] 


[22.80,35.30] 


B{Bs K*-K+) 




[2.56,23.19] 


[2.39,5.78] 


B{B, K-K*+) 




[1.92,6.72] 


[9.73,20.64] 


B{Bs K*-K*+) 




[3.56, 18.76] 


[11.00,45.40] 


A'SpiBs ^ K-K+) 




[0.25,0.49] 


[0.52,0.79] 


A'SfiBskBs ^K*-K+) 




[-0.34,0.07] 


[-0.09,0.64] 


A'SpiB.kBs ^ K-K*+) 




[-0.44,0.05] 


[-0.14,0.61] 


Acf{B, ^ K*-K*+) 




[0.70,0.95] 


[0.79,0.98] 


Ai^'piBs K-K+) 




[0.00,0.06] 


[0.00,0.07] 


A'}}^piB, K*-K+) 




[-0.08,0.02] 


[-0.17,0.00] 


A^SpiBs^ K-K*+) 




[-0.10,0.10] 


[-0.04,0.05] 


A'}^'p{BskBs K*~K+) 




[-0.77,0.27] 


[0.05,0.76] 


A'^Sp{BskBs K-K*+) 




[-0.24,0.76] 


[-0.76,-0.09] 


Acf'iBs K*-K*+) 




[-0.13,0.21] 


[-0.04,0.10] 


fUBs K*-K*+) 




[0.36,0.88] 


[0.75,0.96] 


f^{B, K*-K*+) 




[0.06,0.33] 


[0.02,0.13] 
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Now we turn to the gluino-mediated SUSY contributions to Bg — )■ K^*^ }{(*)+ decays and 
— B^ mixing in the framework of the MI approximation. The following experimental data 



will be used to constrain relevant MI couplings 



io| 



e [0.20,2.84] (at 95% C.L.) 
AM, = 17.77 ±0.12, 
Ar, = 0.19 ±0.07, 



A 



SL = (1.46 ± 0.75) X 10"^ 
= (26.5 ± 4.4) X 10"^ 



(26) 
(27) 
(28) 
(29) 
(30) 



In addition, the same set of the MI parameters also contribute to i? — t- Xg^, which the gluino- 



mediated contribution can be found in Ref. 



28| . Since the experimental measurement of B{B — )■ 



Xs'~f) is in good agreement with the SM expectation, this implies very stringent constraints on 



NP models. We will also use 



a 



B{B X,j) = (3.55 ± 0.24 ± 0.09) x 10" 



(31) 



to constrain the relevant MI parameters. Noted that above experimental data at 95% C.L. will 
be used to constrain the MI parameters. 



3.1 LL insertion 

Let us first consider the LL insertion. The effects of the LL insertions in Bg — )> K^*'^~ K^*^~^ 
decays are almost negligible because there is no the gluino mass enhancement, and B{Bs — >• 
K~K^) given in Eq. fl30l) can not provide any useful constraint on (5^^)23- The bound from 
A^gL is weaker than one from <^f'''^, therefore A^g^ ^i-lso does not give any useful constraint 
when we consider all experimental data given in Eqs. (1261131]) to constrain four kinds of the MI 
parameters. So we only impose the experimental constraints of 5° — 5° mixing and B — t- Xs'^ 
decay, which are shown in Eqs. fl26H28p and Eq. (l?!!) . respectively, to restrict (5^2.)23- 

The constrained spaces of (5^^)23 for rriq = 500 GeV and different x values are demonstrated 
in Fig. m where the allowed parameter space for the MI is shown as dictated by the constraints 
imposed by B{B — ^,7) (yellow), AF, (light gray), cp-^/^'^ (olive) and AM, (pink). The wine 
region shows the allowed regions under the combined constraints of B{B — )■ X,7), B{Bs — )• 
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1 M 




Bound from B{B->Xj) 
Bound from at 
Wf!WWWl| I Bound from (^/** 
i^.^^^^ — Bound from AlVI_. 
iimjiima^i, O Bound from all data 

mm 



0.2 0.4 0.6 

|{8;^),3|forx=0.25 



0.2 0.4 0.6 

|(5^„V I forx-1 
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Figure 4: The allowed parameter spaces of the LL MI parameter constrained by Bs — )■ K~K~^, 
B — and 5^ — 5° mixing at 95% CL. for the squark mass = 500 GeV and the different 
values of x, and denotes the mixing parameters weak phase. 



K~K^), AFs, AMs, A%;i^ and (py^'f'. From Fig. |H we see that the constrained regions are very 
sensitive to the values of x. For x = 0.25,1, as shown in Fig. Il](a-b), the common allowed 
regions are constrained by AF^, (p'^^^'^ and AM^, nevertheless B{B — Xs^) does not give any 
further constraint. For x = 4, 9, we don't show the constraints from B{B Xg'y) in Fig. 
m^c-d) since the whole region of 1(^^^)231 < 1 is allowed by the constraint of B{B — ^^7). As 
displayed in Fig. Il](c-d), the common allowed regions for x = 4 and 9 cases are constrained 
by (j)'^^^'^ and AM^, while AF^ does not give any further constraint. It is worth noting that, 
for X = 0.25, 1,4,9, the lower limit of |(5^2,)23| is also constrained by (p'^^'^'f' since its data are 
not consistent with its SM value at 95% CL. The relevant numerical bounds on |(5^j^)23| with 
different x values are summarized in Table [2l 

49|, the constraint 1(^^^)231 < 0.5 for mg,mq < 600 GeV are derived from B{B — )■ 



In Ref. 



Xg'-f) and B{B — > X^i'^i"). Compared with the existed bound in 49|], for x = 0.25, 1, our upper 
limits of |(5^2.)23| cire at the same order as the previous ones, while the lower limits of 1(^^^)231 
are also given by 0^/'^'''' at 95% C.L.. However, for x = 4, 9, our bounds on 1(5^^)231 are greater 
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Table 2: Bounds on the LL MI parameters from the measurements of B{B — t- Xs7), B{Bs — )■ 
K-K+), ATs, AMs, A'sL and at 95% C.L. for the squark mass rUg = 500 GeV. 



X 


0.25 


1 


4 


9 


MM 


[0.10,0.35] 


[0.22,0.76] 


[0.54,1.00] 


[0.49,1.00] 


0LL (deg.) 


[29,79] 
[-154,-102] 


[24,76] 
[-162,-101] 


[112,170] 
[-69,-10] 


[109,167] 
[-68,-12] 



than ones in Ref. 



1 small tan/3 by AMg, 



50|, for nig = 500 GeV 



49| . Moreover, the bounds on the LL insertion wit 
i3(-B — Xs'-f), A'^cP"' and are also analyzed in detail in Ref. 
and x = 1 case, 1(51^)231 lies in [0.42,0.44] U [0.90,0.95] with tan/3 = 3 and lies in [0.40,0.65] 
with tan/3 = 10. 

The constrained LL MI shown in Fig. H] allows that the theoretical prediction of AMg lies 
in its 95% C.L. experimental range [17.53,18.01] for x = 0.25,1,4,9. However, the ranges of 
AT, and are narrower than their 95% C.L. experimental ranges. For x = 0.25, 1, 
the constrained LL insertion coupling allows ^f/*^"^ G [0.16, 1.26], AF^ G [0.05,0.20] and G 
[0.10,1.00]. For X = 4,9, this coupling allows (pi^^^ G [0.16,0.52], AF, G [0.10,0.20] and 
AI^G [0.10,0.47] . 

Furthermore, we also explore the LL insertion effects in Bg — decays. After 
satisfying all experimental data at 95% C.L. given in Eqs. f l26ll3T]) . the constrained LL insertion 
will not provide significant contribution to B, — )■ if decays. We find the upper limits 

oiB{Bs -» K*-K*+), A'S'piBs^Bs K*-K+,K-K*+), Afp{B, K*-K+) and f^{B, ^ 
K*~ K*^) are slightly decreased from their SM ranges by the constrained LL insertion. The 
lower limits of A^p^iBs — )• K*~ K*'^) and fhiBg K*~K*^) are slightly increased from their 
SM ranges by the constrained LL insertion. The allowed range of Aq^^^^B, — t- K*~ K*~^) is 
increased from its SM prediction [0.70,0.95] to [0.74,1.00] for x = 0.25, [0.77,0.99] for x = 1, 
[0.81,0.97] for X = 4 and [0.76,0.97] for a; = 9, respectively, by the constrained LL insertion. 
While all observables of B, — ?■ if(*)~if(*)+ decays are insensitive to the modulus and weak 
phase of (5^^)23- 
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Figure 5: The allowed parameter spaces of the RR MI parameters constrained by Bs — )■ K~K^ 
decay and 5° — 5° mixing at 95% CL. for the squark mass = 500 GeV and the different 
values of x. 

3.2 RR insertion 

For B — )■ decay, the situation of the RR insertion is very different from the LL one since the 
related NP amplitude (arising from right-handed currents) does not interfere with the SM one. 
Moreover, the effects of the RR insertion in Bg — )■ K~ are almost negligible also because of 
lacking the gluino mass enhancement in the decay. Therefore (5^/^)23 is strongly constrained 
by B^g — Bg mixing. The constrained spaces of {(J'}if>)23 by B^ — B^ mixing for rriq = 500 GeV 
and different x values are demonstrated in Fig. \5[ and the corresponding numerical ranges are 
summarized in Table [31 From Fig. [S] and Table [31 we can see that the allowed moduli and the 
allowed phase ranges of the RR parameters are also very sensitive to the values of x. 

The bound of (5^/j)23 has been obtained in Refs. (49|, |50|. The contributions of the product 
('^LL)23(^^i?)23 are also considered in Ref. j49|, and they obtain |(5^/j)23| < 0.8 for nig, nig < 600 
GeV from B{B X,7) and B{B X,£+r). In Ref. IsQl, the bounds on the RR insertion 
with small tan/3 from AMg, (pi^^'^, B{B — Xg'^), ^cp"* and are also analyzed in detail, 
for ntq = 500 GeV and x = 1 case, |(5^/{)23| lies in [0.36, 0.69] when tan/3 = 3, and there is no 
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Table 3: Bounds on the RR MI parameters from the measurements of Bg ^ K K^, B — )■ 
and B^ — B^ mixing at 95% CL. for the squark mass rriq = 500 GeV. 



X 


0.25 


1 


4 


9 


lfe)23| 


[0.10,0.34] 


[0.23,0.73] 


[0.52, 1.00] 


[0.50,1.00] 


(pRBideg.) 


[20,86] 

[-160,-104] 


[25,75] 
[-153,-101] 


[111,170] 
[-71,-11] 


[118,170] 
[-69,-13] 



common range when tan/3 = 10. 

The constrained RR insertion has the similar effects as the LL insertion on the observables 
of Bg — 7- K^*^~ K^*^^ decays and B'^ — B^ mixing, and we will not show them here. 



3.3 LR insertion 

The effect of the LR insertion is very different from that of either LL or RR. In these decays, 
the LR MI only generates (chromo) magnetic operators Qi'^^sg and Qi^y^g- Especially, the LR 
insertion is more strongly constrained, since their contributions are enhanced by / mi, due to 
the chirality flip from the gluino in the loop. Thus, even a small {^lrIiz can have large effects 
in 5 -> and B, K^*y K'-*^+ decays. 

The constrained spaces of (5^/2)23 from B{B — )■ Xg'j), B{Bs — )■ K~ K^) and 5° — 5° mixing 
for TTiq = 500 GeV as well as different x are demonstrated in Fig. [61 AF^ cannot provide 
any further constraint on (5^/j)23 and we will not show them in Fig. O From the figure, we 
can see that the allowed modulus of the LR MI parameter is very sensitive to the values of x, 
nevertheless the allowed phase range of the LR MI parameter is not changed much for different 
X. We find that B{B — > Xg'j) puts very strong constraints on the upper limit of | (5^^)23 1- And 
0^/"^ also puts very strong constraints on the lower limit of | (5^/^)23! well as the phase of 
(5^^)23. For X = 0.25, 1,4, the allowed spaces from B{B X,7), B{Bs K-K+) and AM, are 
excluded by the constraint from 0^/"^. For x = 9, there is small allowed space from B{B — Xg'y)^ 
B{Bs K-K+), AM, and (pi^^, and it is 1^)23! G [0.08,0.12] U (pLR e [20°, 51°]. 

Previous bound | (5^/2)23! < 0.012 for rrig^mq < 600 GeV has been obtained from the con- 



straint of B{B — Xg'y) in Ref. [49j. Comparing with Ref. [49|], we can see that, as shown in 



Fig. E] (a-c), the bounds for the cases of x = 0.25, 1, 4 from B{B Xg'j) and B{Bs -> K~K+) 
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are stronger than the ones only from B{B — )■ Xg'-f) ahhough they are at the same order. While, 
as shown in Fig. [6] (d) for x = 9 case, the constraint from B{B ^sl) is very strong and 
B{Bs — !■ K~K^) does not give any further constraint. 

Next, we will explore the MI SUSY effects on other observables, which have not been (well) 
measured yet in Bs — )■ decays and B^ — 5° mixing, by using the constrained 

parameter spaces of the LR for x = 9 case as shown in Fig. [6] (d). The numerical results for 
Bs K^*)-K^*)+ and 5° - 5° mixing are summarized in the third column of Table [T] For 
X = 9, the following comments are in order: 

• The LR MI can great increase 

^j/i, f^Qj^ ^Yie SM prediction range [0.034, 0.038] to the 
SUSY prediction range [0.16,0.52], which is however near to the lower limit of the 95% 
CL. measurement. The LR MI has been restricted by the experimental upper limit of 
B{Bs — !■ K^K^), and the allowed range of B{Bs — > K^K'^) is significantly shrunken 
from its SM prediction [9.20,45.52] x 10^^ to [22.80,35.30] x 10^^ by the constrained LR 
insertion. 




[{5' J J for lis^JJIorx^d 
(c) (d) 

Figure 6: The allowed parameter spaces of the LR MI parameters constrained by Bg K~K^ 
decay and 5° — -B^ mixing at 95% CL. for the squark mass = 500 GeV and for the different 
values of x. 
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• The constrained LR could affect the branching ratios significantly. The allowed upper 
limit of B{Bs — ?■ K*~ K^) could be reduced from its SM prediction, and the allowed 
values of B{Bs — t- K~K*~^, K*^K*^) are great increased by the constrained LR insertion. 
The range of SUSY prediction of B{Bs — )■ K~ K*^) could differ from its SM expection 
significantly. 

• The constrained LR insertion has great contributions to all mixing CPAs in Eg — )■ 
Ki*)-K^*)+ decays, and all mixing CPAs could be largely enhanced. In addition, the 
constrained LR insertion could change Acp{BshBs K*~K~^,K~K*~^) a lot. 

• The polarization fraction fiiBs -> K*~K*~^) can be enhanced much by the constrained 
LR insertion. 

For LR insertion with x = 9, we can present the distributions and correlations of B, A(^<p, 
Aq'^p, fL,± within the modulus or weak phase of the constrained LR MI parameter space in Fig. 

by two-dimensional scatter plots. The LR MI effects on all observables of Bg — > 
decays are displayed in Figs. (I7J|8]). Fig. [7] and Fig. |8]show the sensitivities of all observables 
to 1(5^^2)23! and (pLR, respectively. In addition, for comparing conveniently, we show the SM 
bounds of these observables by orange horizontal dash lines and the limits of the measurements 
of B{Bs — 7- K~ K^) at 95% C.L. by the cyan horizontal solid lines. From Fig. [7] (a-d) and Fig. 
[8] (a-d), one can find that B{Bs — )■ K*~ , K~ K*~^) have mild sensitivities to both 1(5^^)231 
and (pLR, while B{Bs — )■ K*~K*~^) is insensitive to |(5l/j)23| or (p^p. As shown in Fig. [7](e-h) 
and Fig. [He-h), the LR insertion has positive effects on all four mixing CPAs, and they are 
sensitive to both | (5^^)23! 'Plr- So the future measurement of any mixing CPA could further 
restrict both |(5i^)23| and Fig. [8](i) and (k) show ^cp(^s ~^ K-K+,K-K*+) are mildly 
sensitive to (f)LR. Fig. [7] (m-n) and Fig. E (m-n) display that A'^Sp^BgkBs -> K*^K+, K~K*+) 
are sensitive to both \{S'lp)23\ and (pLR- As for the LR insertion effects on — ^ K*~K*~^) 
and f±{Bs — )■ K*~K*~^), we show them in Fig. [7] (o-p) and Fig. [8] (o-p), we can see 
K*-K*+) and f\{Bs K*- K*+) could be affected significantly by the LR MI. 
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3.4 RL insertion 

The SUSY contributions of the RL insertion also pick up an mg/mb enhancement relative to 
the SM. Compared to the LR case, the RL situation is very different since the related NP 
amplitude does not interfere with the SM one in B{B — )■ Xs-f). The RL insertion is much more 
strongly constrained by B{B — > Xs'~f). 

The constrained spaces of (5^^)23 from B{B — )■ X^'j), B{Bs — K~ K^) and 5^ — mixing 
for TTiq = 500 GeV and different x are demonstrated in Fig. O AF^ can not provide any further 
constraint on (5^x)23 which is not shown in Fig. [9l As shown in this figure, there is no common 
space from B{B — > Xg'y), B{Bs K~K^), AMg and 0^/"^ since B{B Xs'y) puts very strong 
constraints on the upper limits of 1(4^)231, roughly |(5^ij)23| < 0.0057,0.0086,0.020,0.036 for 
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Figure 7: The effects of (5^/2)23 for x = 9 case in Bg — ^ decays. B are in units of 10~^. 

The orange horizontal dash-dot lines denote the limits of SM predictions, and the cyan horizontal 
solid lines represent the 2a error bar of the measurements. (The same in Fig. [8|). 



20 





50 




40 






A 


30 


1^ 

cd" 


20 


CD 






10 




Figure 8: The effects of cpLR ^oi x = 9 case in Bs K^*)-K'^*^+ decays. 



X = 0.25,1,4,9, respectively. So the RL insertion cannot accommodate the current data of 
B{B X,7), B{Bs K-K+), AM, and (pi^^ simultaneously. 



4 Conclusions 

Motivated by the recent measurements from CDF and D0 collaborations, we have studied the 
gluino- mediated SUSY contributions to 5°— 5° mixing, B, and B X^'j decays 

with the MI approximation. Considering the theoretical uncertainties and the experimental 
error bars, we have obtained fairly constrained parameter spaces of LL, RR, LR and RL Mis 
from the present experimental data of B^ — B^ mixing, Bg — t- K~K^ and B — )■ Xg'^ decays. 
Furthermore, using the constrained MI parameter spaces, we have predicted the MI SUSY 
effects on the observables of four Bg — )■ K'^*^' K^*'^^ decays, which have not been measured yet. 
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For the LL and RR Mis, the strong constraint arises from — mixing, and B{Bs — >• 
K~K'^) as well as B{B — )■ Xg'y) cannot provide any further constraint on {^ll Riil'^'i- 
have found that, for x = 0.25, 1,4,9 cases, the constrained LL and RR Mis have little effect 
on the observables of Bg K^*)- K^*)+ decays. The upper limits of B{Bs K*~K*'^), 
A^'^iBskBs K*-K+,K^K*+), ^Sp(^« ^ K*^K+) and f±{Bs K*~K*+) are slightly 
decreased from their SM values. The lower limits of Aqp^{Bs — )■ K*^K*'^) and 
K*~K*~^) are slightly increased from their SM values. The allowed range of ^^p*"(i?s — )• 
K*~K*~^) is enlarged. 

For the LR and RL Mis, B{B — )■ Xg'y) puts particularly strong constraints on the upper 
limits of |(^L/j.ijx)23|, and 0^/"^ also puts very strong constraints on the lower limits of KSlb. rl)23| 
as well as the phases of (5^^ rl)23- So only very narrow space of the LR MI for x = 9 case could 
explain the 95% C.L. experimental data of AF,, AM,, 0^/'^'^, B{Bs K- K^) and B{B ^ 
Xs7) simultaneously. We have found the constrained LR insertion for x = 9 still have sizable 
effects on all observables of Bs — )■ decays except Acp{Bs — )■ K~K^). In addition, 

we have presented the sensitivities of the observables to the constrained LR parameter spaces in 




Bound from B(B->X j) 

1 Bound from B(B^->K^K') 
- Bound from AM 

s 

Bound from 




1(6" for x=4 K5V),3lforx=9 
(c) (d) 

Figure 9: The allowed parameter spaces of the RL MI parameters constrained by Bg — )■ K~K^ 
decay and 5° — i?, mixing at 95% C.L. for the squark mass rriq = 500 GeV and the different 
values of x. 
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Figs. OlHi We have found that all mixing CPAs of Bs K^*)-K^*)+ are very sensitive to both 
|((5i^)23| and moreover, K*~ K+ , K*+), A^^^Bs ^ K*-K+), A^^UBs^Bs ^ 

K*~ , K' K*'^) and fL,±iBs K*~ K*^) have some sensitivities to (^^^2)231 or (piR- So the 
future measurement of any mixing CPA could be very useful to shrink/reveal/rule out the 
relevant LR MI parameter space. The results could be useful for probing SUSY effects and 
searching direct SUSY signals at Tevatron and LHC in the near future. 
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Appendix: Input parameters 

The input parameters are collected in Table |H We have several remarks on the input parame- 
ters: 

• Wilson coefficients : The SM Wilson coefficients Cf are obtained from the expressions 



in Ref. 



24|. 



CKM matrix element : For the SM predictions, we use the CKM matrix elements from 
the Wolfenstein parameters of the latest analysis within the SM in Ref. 52|, and for 



the SUSY predictions, we take the CKM matrix elements in terms of the Wolfenstein 



parameters of the NP generalized analysis results in Ref. 52 1. 



Masses of SUSY particles: When we study the SUSY effects, we will consider each possi- 
ble MI (5^^)23 for AB = LL, LR, RL, RR only one at a time, neglecting the interferences 
between different insertions products, but keeping their interferences with the SM am- 
plitude. We fix the common squark masses rrig = 500 GeV and consider four values of 
X = 0.25, 1, 4, 9 (i.e. m-g = 250, 500, 1000, 1500 GeV) in all cases. 
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Table 4: Values of the theoretical input parameters. To be conservative, we use all theoretical 
input parameters at 68% CL. in our numerical results. 



mw = 80.398 ± 0.025 GeV, m^^ = 5.366 GeV, m^_^^ = 0.892 GeV, m^^ = 0.494 GeV, 




mt = 171.3ii;6 GcV, mfc(mb) ^ (4.20 ±0.07) GcV, m^(2GcV) = (0.105lJ^;;^3^) GcV, 




r,^ ^ (1.472l0;0^26) 


[5JJ 


The Wolfenstein parameters for the SM predictions: 




A = 0.810 ± 0.013, A = 0.2259 ± 0.0016, p = 0.154 ± 0.022, 77 = 0.342 ± 0.014. 




The Wolfcnstcin parameters for the SUSY predictions: 




A = 0.810 ± 0.013, A = 0.2259 ± 0.0016, p = 0.177 ± 0.044, 77 = 0.360 ± 0.031. 


[52] 


]k = 0.160 GeV, /a- = (0.217 ± 0.005) GeV, = (0.156 ± 0.010) GeV, 




jsB^^K- (0) = 0.360 ± 0.034, Af (0) = 0.233 ± 0.022, A^'^^' (0) 0.181 ± 0.025, 




y^=^-^*(0) =0.311 ± 0.026, F^^^^iS>) = 0.30t°;°^. 


[53, 5^ 


= (0.245 ± 0.025) GeV, /s^ ^ Bb, = 0.270 ± 0.030 GcV. 


[55] 


?72B = 0.55 ±0.01. 


m 


af = 0.2 ± 0.2, af = 0.1 ± 0.3, af * = 0.06 ± 0.06, af * = 0.1 ± 0.2. 


05,26] 


S(^)(mb) = 0.86(2) (It), S(^)(mfc) =0.83(2)(4), ^^^^(TTZfa) = 1.03(4)(9), 




B'^:\m,) = 1.17(2) {X^) > - 1-94(3) (tf ). 


[57] 
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